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ABSTRACT 



A waterflood of the Second Wall Creek Sand of Naval Petroleum 
Reserve No. 3 is being considered in order to increase the ultimate oil 
production. A streamline model of the Second Wall Creek Sand was 
developed here and used to predict the response of the reservoir to a 
waterflood. 

Results of the model study indicate that the oil production of the 
Second Wall Creek Sand can be increased by as much as 15, 000, 000 
barrels over a twelve year period by waterflooding. This would bring 
the total production available from the Second Wall Creek Sand to around 
32, 000, 000 barrels. 
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CHAPTER I 



INTRODUCTION 

Naval Petroleum Reserve No. 3 is located in Natrona County, 
about thirty-five miles north of Casper, Wyoming. The major geo- 
logical structure on the reserve is an anticline called "Teapot Dome". 
The Second Wall Creek Sand is the most important of the several oil 
bearing formations on the anticline. 

The Second Wall Creek Sand is a solution gas drive reservoir. 
There appears to be very little natural water drive. Early estimates, 
based on recovering one -seventh of the original oil in place, predicted 
the recoverable reserves from this formation to be around 17,000,000 
barrels (1). 

Waterfloods of the same formation in neighboring fields have 
proven to be successful in increasing production. Therefore, in order 
to enhance the recovery, the development and production of the Second 
Wall Creek Sand as a waterflood is being considered. 

The objective of this work is to apply streamline simulator 
techniques to the Second Wall Creek Sand and predict the results of a 
waterflood under certain operating conditions. 
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CHAPTER II 



THE STREAMLINE MODEL 

The theory and development of the streamline model have been 
rigorously described by LeBlanc and Caudle(2), Rust(3), Wessels{4) 
and Kazmann(5). A review of the assumptions and the mathematics 
is presented here. 

2. 1 Basic Assumptions 

The basic assumptions used in formulating the streamline model 
are as follows : 

1. Negligible gravitational effects within the reservoir. There 
is no gravity segregation of fluids. 

2. Steady-state flow within the reservoir. Incompressible 
reservoir fluids. Operations are at nearly constant pressure. 

3. Homogeneous, isotropic reservoir. 

4. Uniform dispersion of free gas. 

5. Oilbank formation. 

The simulator used is designed specifically for oil bank build-up 
and will not accurately predict the results of other operations such as 
pressure maintenance. 
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2. 2 Fundamental Equations 

The starting point for the development of the streamline model is 
the continuity equation. This partial differential equation is derived from 
a three dimensional material balance and is written as: 



S(^P 



» ) 



dx 




( 2 . 2 . 1 ) 



The volumetric flux may be expressed by Darcy's Law in the 
differential form as: 



q_ _k_ d $ 

A " " \s dl 



( 2 . 2 . 2 ) 



When the volumetric flux components in the x, y and z directions, as 
described by Darcy's Law, Equation 2.2.2, are substituted into 
Equation 2.2. 1, the following equation results: 

a / pk x a$\ a / pk y a$\ , _a_ f pk z a$\ , bp 

bx \ |i bx / by ' (J. by) bz \ |-l bz) bt (2.2.3) 

By using the assumptions listed previously. Equation 2.2.3 
may be simplified. The resulting equation may be written in the form 
of a La Place equation. 

a 2 p , * 2 p _ 
a * 2 a y 2 
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The solution of Equation 2.2.4 yields the pressure at any point 
in a porous medium for which the assumptions made are valid. Knowing 
the pressure at all points in a porous medium, the average velocity of a 
fluid particle, and thus its movement, may be computed. 

By transforming to cylindrical coordinates and allowing radial 
flow about a point, Equation 2. 2.4 may be solved for the flow potential 
in a horizontal, isotropic, homogeneous porous medium. Restoring 
the resulting equation to rectangular coordinates and using the rule of 
superposition for a multi -well system, the solution of Equation 2.2.4 
becomes 



p(x,y) = p 
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(2.2. 5) 



where p is the mean reservoir pressure, 
m 

In a porous medium, the relationship of flux to velocity is 



v = 



u 

? 



k dp 
4>li dl 



( 2 . 2 . 6 ) 



By differentiating Equation 2. 2. 5 with respect to x and substituting the 
resulting derivative into Equation 2.2. 6 the velocity in the x direction 
in a rectangular coordinate system becomes: 



1 n x - x 
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In a similar manner, the velocity in the y direction becomes: 



v y (x, y) = 



2rrcf>h 



n 

E 



y - y i 



i = 1 * ( x - x i) 2 + ( y - y i )‘ 



( 2 . 2 . 8 ) 



Using vector addition, the velocity of a particle of fluid is: 



2 2 
V = V + V 

vj x y 



(2.2.9) 



2. 3 Finite Difference Equations for Movement 

In order to apply Equations 2. 2. 7 and 2. 2. 8 to the streamline 
model, it is assumed that the velocity remains constant during any one 
small distance increment, As. Then the time increment. At, is 
calculated by: 

A s 

At = — (2.3.1) 

v 

From the particles starting point (x , y ), it will move to a new point 

s s 

calculated as follows : 



x = X + V (x ,y )At 

s + 1 s x s s 



y . 1 = y + v (x ,y )At 
s+1 s v ss 



(2.3.2) 

►At 

y ' s ’ ' s 

By starting on the wellbore of an injection well and successively 
applying Equations 2.2.7 through 2.3.2, a trace of any particle's move- 
ment may be made. This trace, called a streamline, is followed until it 
reaches a production well whereupon it is terminated. The travel time 



for a single streamline is the summation of the time increments 
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corresponding to the number of steps taken along the streamline. 

2. 4 Streamtube Concept 

In actual computation, a representative number of streamlines are 
chosen to represent all fluids emanating from the source. An example of 
the streamlines for a bounded pattern are shown in Figure 1. A few 
typical streamlines are shown in Figure 2. By dividing the distance 
between two adjacent streamlines into two equal parts, the streamline 
may be said to represent the flow of all the particles within the dividing 
lines. A streamtube has thus been described and is shown in Figure 3(3). 
As a fluid particle is traced along its streamline, the volume of fluid 
produced at the streamline’s production well can be calculated and 
accumulated. 

2. 5 Multiple Fluid Flow 

The streamline model described so far deals with a single fluid. 

In a waterflood operation, water is injected into the reservoir to displace 
the oil and gas. Therefore as many as three fluid regions may be present 
in the reservoir at any time and may each have a different mobility. 

The first assumption for the multiple fluid system is that the oil 
is miscibly displaced by the water in a piston-like manner. The oil is 
banked up ahead of the water such that all oil displacement takes place at 
the leading edge of the water front. A diagram of the concentrations is 
shown on page 9. 
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Figure 1 Streamlines for a Bounded Pattern. 
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Figure 2 Typical Streamlines. 
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Figure 3 Typical Streamtube. 
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With the assumption of piston-like displacement the displacing fluid inter- 
face can be easily followed along the center pathline of a streamtube. 

The most significant approximation used in extending the stream- 
line model for multiple fluids is that the pathlines will be the same as the 
streamlines generated for a single fluid. This approximation, while not 
correct, has been found to be reasonable except at extreme values of 
mobility ratio (3). 

As the fluids proceed along the streamtube the flow rate of the 
fluid stream will change because of the different fluid region mobilities. 
To correct for the changing fluid flow rate the concept of conductivity 
ratio is applied . Conductivity and conductivity ratio are developed and 
explained by Caudle (6). Briefly the conductivity ratio is a comparison 
of the conductivity of the medium at some time when more than one fluid 
is flowing to the conductivity of the medium initially when only one fluid 
was present. The general equation for the conductivity ratio is 
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(2.5. 1) 



The velocity Equations 2.2.7 and 2.2.8 are multiplied by the conduc- 
tivity ratio to correct for multiphase flow. 



2. 6 Modification for Thickness Variation 

An empirical modification to allow for varying thickness has been 
developed by Wessels(4). This modification is useful if the thickness at 
each point in the reservoir can be defined. A system of grid squares 
with assigned thicknesses is used to describe the reservoir thickness 
variation. 

The modifications for thickness variation are as follows: 
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(2.6.3) 



The factor (h^/h.) is an approximation factor which reflects an increasing 
velocity as the thickness decreases and a decreasing velocity as the 
thickness increases. The diagram below shows the relationship between 
reservoir thickness and velocity for a single streamline. 
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2. 7 Mathematical Bounding of the Streamline Model 

An imaging technique developed by Lin (7) which can determine 
the necessary rates for a preselected set of image wells that will 
confine a system of streamlines within a given area was used for this 
study. The technique uses a least squares method to determine the 
image well flow rates so that a no -flow boundary is established between 
each of the reservoir bounding points. 



CHAPTER III 



APPLICATION OF THE STREAMLINE MODEL TO THE 
SECOND WALL CREEK SAND 

3c 1 Reservoir History 

President Wilson's Executive Order of April 30, 1915 designated 
the Teapot Dome area in Wyoming as Naval Petroleum Reserve No. 3 
(NPR 3)o On April 7, 1922, the reserve was leased to the Mammoth Oil 
Company for the purpose of exploitation* Following a decision of the U.S. 
Supreme Court all of the producing wells on the reserve were shut-in on 
December 31, 1927. The field remained shut-in until 1951-1953 when an 
exploratory program was initiated; thereafter it was shut-in until 1958 
when an offset drilling program was instituted to protect against drainage 
by adjacent operators. 

The wells are distributed unevenly in the Second Wall Creek Sand 
because development was started at the north end of the reserve and was 
stopped prior to completion. Based on approximate surface area the 
northern third of the reserve has 69 wells, the middle third has 23 wells 
and the southern third has 13 wells. 

Approximately 3, 600, 000 barrels of oil were produced prior to 
December 31, 1927. 1 , 200, 000 barrels of oil have been produced since 

1951 with most of this production occuring after 1958. 
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3.2 Location of Oil and Gas 

Figure 4 shows the approximate locations of the gas cap, oil zone 
and water in the Second Wall Creek Sand (8) . 

The Second Wall Creek Sand just beyond the northwestern boundary 
of NPR 3 has been under waterflood for several years, The Navy has been 
producing wells along this boundary to protect against drainage. Conse- 
quently a portion of the Naval Reserve has been swept by the waterflood. 

To allow for this production, 17 0 acres along the northwestern boundary 
were omitted from the area which was modeled. The total area modeled 
was 3000 acres. 

3. 3 Reservoir Properties 

Reservoir data and operating conditions were obtained from 
Tesoro Petroleum Corporation and U. S. Navy drawings of NPR 3 
(8, 10, 11), Geological Survey and U. S. Navy publications (1, 9), 

U. S. Navy correspondence (12, 13) and various core data and logs pro- 
vided by the Officer in Charge of Naval Petroleum Reserve No. 3. 

Table I lists the reservoir and fluid properties used for the 
modeling study. A discussion of significant individual items follows. 

Porosity and Saturations The average porosity and connate water 
saturation used are the same as those used in previous estimates of 
waterflooding at the north end of the reservoir (12, 13). No significant 
difference was noted between these figures and the core data examined. 



* 
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Figure 4 Location of Oil and Gas in the Second Wall Creek Sand. 



TABLE I 



RESERVOIR AND FLUID PROPERTIES 

Porosity (Avg), % 17.5 

Connate Water Saturation, % 35 

Gas Saturation, Initial, % 5 

Gas Saturation, Residual, % 0 

Oil Saturation, Initial, % 60 

Oil Saturation, Residual, % 20 

Major Fault Blocks 10 

Net Thickness, Ft. 21-60 

Reservoir Area, Acres 3000 

Depth (Avg), Ft. 2800 

Temperature, °F 140 

Pressure, PSI 1100 

Permeability (Avg). MD 10 

Formation Volume Factor, Oil 1.2 

Formation Volume Factor, Water 1.0 

Oil Gravity, °API 3 8 

Oil Viscosity, CP 6 

Water Viscosity, CP . 7 

Gas Viscosity, CP . 12 

Water/Oil Mobility Ratio 4 

Gas /Oil Mobility Ratio 65 

Gas Oil Ratio, SCF/BBL 1000 

Solution Gas Oil Ratio, SCF/BBL 50 
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For this study it was assumed that the waterflood would start at 
the very beginning of the production. The model was run with an initial 
gas saturation of 5% to allow oil bank build-up to occur and to show the 
effect of fill-up time. The initial oil saturation used was 60%. 

The residual gas saturation was set at zero because it was 
assumed that all of the free gas would either be displaced or put back 
into solution as the oil bank passed. 

The residual oil saturation after the waterflood passed was set 
at 20%. This is almost twice the 12% average residual oil saturation 
found in the cores available. The residual oil saturation in the cores 
was considered an absolute minimum which would be difficult to achieve 
in the reservoir. 

Fault Blocks The Second Wall Creek Sand is heavily faulted and 
was, for the presentation of well productivity data, divided into about 
twenty possible fault blocks (12). In order to model the field, taking into 
account the faulting, the field was divided into ten sections along the 

t 

largest faults. It was assumed that the faults were sealed and the 
sections were then modeled independently. This had the added advantage 
of making the modeled reservoirs smaller and easier to manage. Also, 
properties which vary within the reservoir could be easily changed from 
section to section. Figure 5 shows the configuration of the ten fault 



blocks chosen. 
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Figure 5 Fault Blocks Used in Modeling the Second Wall Creek Sand. 
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Net Sand Thickness Thin beds of shale are common throughout 
the Second Wall Creek Sand, but they thicken toward the south and east 
and split the sand into several distinct layers. Some wells have three or 
more layers of sand separated by shale beds from two to twenty-five 
feet thick(12)o 

Combining the information available from the cores with the 
general trends mentioned above, a percentage discount factor was used 
for each fault block to reduce the gross sand thickness to a net sand 
thickness. Table II shows the discount factors that were applied. 

Figure 6 shows the gross thickness of the Second Wall Creek Sand(10). 

Permeability The temperature, pressure and permeability 
were used to determine the fluid and fluid flow properties. The perme- 
ability itself is not used in the streamline model calculations because the 
calculations are made at a constant flow rate. This means that the 
injection pressure will have to be great enough to make the injectors 
flow at the required rate. Rate selection will be discussed in the 
next chapter. 

Gas -Oil Ratios Gas -oil ratios vary greatly depending on where 
in the field the well is located. The gas -oil ratio and solution gas -oil 
ratio are used by the computer program to indicate the arrival of the 
oil bank. Representative values of the two ratios were used but they do 
not represent a rigid material balance. 
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TABLE II 

DISCOUNT FACTORS APPLIED TO 
GROSS SAND THICKNESS 
TO OBTAIN NET SAND 
THICKNESS 
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Figure 6 Isopachous Map of the Second Wall Creek Sand- 
Gross Sand Thickness. 
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3. 4 Operating Conditions 

Table III is a listing of the operating conditions used for the 
simulation. 



TABLE III 

OPERATING CONDITIONS FOR THE WATERFLOODING 
OF THE SECOND WALL CREEK SAND 



No. of Production Wells 

No. of Injection Wells 

Well Pattern 

Well Spacing, Acres 

Maximum Oil Production Rate, 
STB/Well/Day 

Oil Production Rates During Fill Up, 
STB/Well/Day 

Total Water Injection Rate, STB/Day 



135 

132 

5 -Spot/ Line Drive 
10 

46-300 

3-25 

23,000 



Field Development Based on a proposed Navy Development 
Plan(ll) the field was assumed to be drilled on ten acre spacing. 

Well water injectivity was assumed to be in the same range as 
the well's initial oil productivity. Therefore, approximately a one-to- 
one ratio of injectors to producers was maintained. 

The injection-production well pattern used was a combination 
five -spot and line drive. The injection wells were generally placed 
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along the outside, nearest the oil/water contact, in a line drive pattern to 
force the oil toward the center. Then, where area allowed, injection wells 
were placed within the producers in a five -spot pattern. 

Production Well Rates In the streamline model the oil production 
rate of a well during the fill-up time is set equal to the present rate to 
which the well's oil production has declined. After oil bank breakthrough 
the oil production rate of a well is set equal to the initial or maximum 
oil rate of that well. 

The initial rates for each production well were determined from 
a composite drawing of the Second Wall Creek Sand production(8). Since 
the reservoir will be waterflooded from the start, the oil rates during 
fill-up were set at 10% of each well's initial rate. At reservoir conditions, 
the total water injection rate was set equal to the total oil production rate 
after fill-up to satisfy the steady state assumption. 

The individual injection rates were assigned to the injection wells 
in such a manner to contain the streamlines as much as possible to the 
oil bearing area. Some flow into the gas cap and the watertable cannot 
be stopped but proper balancing minimized the losses. 

Figures 7 and 8 show the well pattern used and the streamlines 
generated by the model reservoir. See Appendix A for a listing of the 
individual well.data used. 
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Figure 7 Model Generated Streamlines -Northern Half of 
Second Wall Creek Sand. 



24 




Figure 8 Model Generated Streamlines-Southern Half of 
Second Wall Creek Sand. 



CHAPTER IV 



RESULTS OF THE STREAMLINE MODEL SIMULATION 

4. 1 Predicted Production History 

Figures 9, 10 and 11 show the predicted response of the Second 
Wall Creek Sand to a waterflood operated as previously described. It is 
predicted that 18, 000, 000 barrels of oil will be recovered in five years 
and 30,000,000 barrels of oil after ten years. Water breakthrough occurs 
around 1, 000 days and the water -oil ratio climbs to over five after ten 
years. Ultimate recovery at a water-oil ratio of twenty is predicted 
to be around 32, 000, 000 barrels and takes twelve years. 

4. 2 Effects of Major Assumptions on Recovery Predictions 

Although the oil recovery predicted agrees with the present 
reserve estimate of 28, 800, 000 barrels (9), the predictions of both 
quantity and time are only as good as the assumptions made. 

In order to provide a proper perspective, the variables and 
assumptions which could most effect the results will be discussed. 

Heterogeneity The model was run assuming that the reservoir 
sand in each fault block was homogeneous. Frequently, however, hori- 
zontally bedded sandstones exhibit different permeabilities among the 
individual stratum. The presence of high permeability strata will 
reduce the efficiency of the waterflood. The high permeability zones 
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Water Injection Rate = 23, 000 STB/DAY 
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Figure 9 Oil and Water Production Rates During Waterflood Second Wall Creek Sand. 
Predictions Made Using Streamline Model. 
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Figure 10 Oil and Water Production During Waterflood Second Wall Creek Sand. 
Predictions Made Using Streamline Model. 
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Figure 11 Water Injected and Total Liquids Produced During Waterflood Second Wall 
Creek Sand. Predictions Made Using Streamline Model. 
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water-out faster than the rest of the reservoir and then increase the 
water oil ratio much sooner than expected. If high permeability strata 
are shown to be a significant factor the predictions presented here will 
be optimistic. 

Initial Gas and Oil Saturation The presence of free gas in the 
oil zone of the reservoir means that there will be a certain time lag, 
fill-up time, before the oil bank will reach the production wells. This 
time lag is approximately equal to the time required to inject enough 
water to fill the gas space and the space of the oil that is produced. The 
initial gas saturation chosen to represent the reservoir was 5%. If the 
quantity of gas present in the reservoir is actually less than 5% the oil 
production rate will increase earlier and if the space is assumed to 
be filled with oil, more oil will be recovered. If the gas saturation is 
greater than 5% the opposite will be true. 

Production and Injection Rates The effect of increasing or 
decreasing the injection and production rates would be to condense or 
expand the time scale. If for example it was found that the reservoir 
could operate at injection and production rates 50% greater than those 
assumed here, the recovery that would have taken 1, 500 days could be 
achieved in 1,000 days. See Appendix B for an example of rate adjust- 



ment to obtain uniform water breakthrough times. 



CHAPTER V 



CONCLUSIONS AND RECOMMENDATIONS 
5. 1 Conclusions 

The waterflood of the Second Wall Creek Sand will increase the 
total oil production to around 32, 000, 000 barrels. This is 15, 000, 000 
barrels over that which could be expected from solution gas drive primary 
production. At the operating rates used the waterflood will require around 
twelve years. 

5. 2 Recommendations 

1. It is recommended that the waterflooding of the Second Wall 
Creek Sand be planned and cost evaluations made. The streamline model 
developed here is well suited for studying development and operating 
considerations so that project economics can be evaluated. 

2. Since the effects of high permeability strata on secondary 
recovery operations are always adverse and sometimes severe, it is 
recommended that the possible existence of such strata be thoroughly 
investigated during waterflood planning and development. 
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NOMENCLATURE 



Symbol 

A 

h 

k 

1 

n 

P 

q 

As 

s 

At 

t 

u 

V 

x, y, z 

y 

x 

u 

p 

o 

$ 



Definition 

Cross Sectional Area 

Reservoir thickness 

Permeability 

Length; number of fluids 

Number of wells 

Pres sure 

Flow rate 

Distance increment along path line 
Saturation 

Time step along path line 

Time 

Flux 

Average fluid particle velocity 
Rectangular coordinates 

Conductivity ratio 

Fluid mobility 

Viscosity 

Density 

Porosity 

Potential 
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Subscripts 

cw 

gi 

g r 

i 

j 

k 

m 

0 

01 
or 
s 

P 

w 

X 

y 

z 



Connate water 

Initial gas 

Residual gas 

Index of wells; in 

Index of particle location 

Index of fluids 

Mean 

Oil; out 

Initial oil 

Residual oil 

Starting point of fluid particle 
Production 
W ater 

In the x-direction 
In the y -direction 
In the z -direction 



APPENDICES 



APPENDIX A 



INDIVIDUAL WELL DATA 
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TABLE IV 

WELL DATA— FAULT BLOCK 1 
Refer to Figure 12 



Well 

No. 


Rate, STB/Day 
(+Inj , -Prod) 


Thickness 

Ft. 


Well 

No. 


Rate, STB /Day 
(+Inj, -Prod) 


1 


87. 


29 


21 


96 


2 


54. 


27 


22 


-96 


3 


-90. 


30 


23 


48 


4 


52. 


26 


24 


-72 


5 


-87. 


29 


25 


-72 


6 


93 


31 


26 


-78 


7 


56 


28 


27 


84 


8 


-90 


30 


28 


-90 


9 


-96 


32 


29 


96 


10 


78 


26 


30 


-96 


11 


-81 


27 


31 


-96 


12 


90 


30 


32 


48 


13 


-96 


32 


33 


-72 


14 


96 


32 


34 


72 


15 


48 


24 


35 


-72 


16 


72 


24 


36 


78 


17 


-78 


26 


37 


-84 


18 


-81 


27 


38 


90 


19 


87 


29 


39 


-96 


20 


-93 


31 


40 


96 



Thicknes s 
Ft. 

32 

32 

24 

24 

24 

26 

28 

30 

32 

32 

32 

24 

24 

24 

24 

26 

28 

30 

32 



32 



TABLE IV (Continued) 



Well 

No. 


Rate, STB/Day 
(+Inj , -Prod) 


Thickness 

Ft. 


Well 

No. 


Rate, STB/Day 
(+Inj, -Prod) 


Thicknes s 
Ft. 


41 


48 


24 


61 


-90 


30 


42 


-72 


24 


62 


81 


27 


43 


72 


24 


63 


48 


24 


44 


-72 


24 


64 


36 


24 


45 


72 


24 


65 


72 


24 


46 


-78 


26 


66 


-78 


26 


47 


84 


28 


67 


93 


31 


48 


-90 


30 


68 


-90 


30 


49 


96 


32 


69 


84. 


28 


50 


-96 


32 


70 


-84 


28 


51 


96 


32 


71 


48. 


24 


52 


48 


24 


72 


48. 


24 


53 


-72 


24 


73 


-78 


26 


54 


-72 


24 


74 


-78 


26 


55 


-72 


24 


75 


36 


24 


56 


90 


30 


76 


40 


24 


57 


-90 


30 


77 


48 


24 


58 


90 


30 








59 


-90 


30 








60 


90 


30 









Production well oil rates during fill-up equal 10 STB/Day 
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\ 



39 28 




/z 


7 


y 


\ 


O X 


O 


X 


X 


X 


90 3 £> 


2/ 


/3 


8 




2 


X © 


y 


o 


O 


o 


X 




22 


/y 


9 


(p 


3 




O 


X 


o 


X 


O 






28 









\ 



\ 



27 



x Injection Well 
O Production Well 



FAULT BLOCK 2 









\ 35" 29 2S ZG 

\ x o x o 

x "o S' 

\ 3 J o' 

\ 38 
N X 


2^ 22 Zo 17 13 8 S' / 

x oxoxoox 

V If 2 o 3 *' o 8 'l '$ % 

V \i’ 7 'i 

a ; x 


\ 

2 \ 

X a \ 

O X \ 

'h % *X \ 


Section 29 


28 


27 



Figure 12 Well Locations for Fault Blocks 1 and 2. 
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TABLE Y 

WELL DATA— FAULT BLOCK 2 



Refer to Figure 12 



Well 

No. 


Rate, STB/Day 
(+Inj, -Prod) 


Thicknes s 
Ft. 


Well 

No. 


Rate, STB /Day 
(+Inj , -Prod) 


Thicknes s 
Ft. 


1 


* 


28 


20 




29. 4 


2 


* 


30. 1 


21 


❖ 


30. 8 


3 


* 


31.5 


22 


-175 


29. 4 


4 


❖ 


32. 2 


23 


-180 


30. 1 


5 


-175 


28.7 


24 




28. 7 


6 


-180 


30. 1 


25 


4 


29. 4 


7 


-195 


32. 2 


26 


-155 


25. 9 


8 


-205. 


34.3 


27 


-175 


28.7 


9 


-175 


28.7 


28 




24. 5 


10 


-185 


30. 8 


29 


-135 


22.4 


11 


-200 


32.9 


30 


-147 


24. 5 


12 


-210 


35. 


31 


* 


25. 9 


13 




29.4 


32 


-160 


26. 6 


14 


❖ 


30. 8 


33 




28. 7 


15 




33. 6 


34 


-175 


29.4 


16 


'i" 


35. 


35 




21. 0 


17 


-175 


29. 4 


36 


❖ 


23. 1 


18 


-185 


30. 8 


37 


sjc 


25. 9 


19 


-200 


32. 9 


38 


* 


28. 7 



Production well oil rates during fill-up equal 20 STB/Day. 

^Program calculates injection well rates based on reservoir thickness 
and steady state material balance. 



TABLE VI 



WELL DATA— FAULT BLOCK 3 



Refer to Figure 13 



Well 

No. 


Rate, STB /Day 
(+Inj, -Prod) 


Thickne 

Ft. 


1 


140 


35. 


2 


140 


35 


3 


140 


35 


4 


-178. 5 


35.7 


5 


-182 


36.4 


6 


-178. 5 


35.7 


7 


218.4 


36.4 


8 


226.2 


37.7 


9 


226. 2 


37. 7 


10 


-182 


36.4 


11 


-192. 5 


38. 5 


12 


-203. 


40. 6 


13 


-188. 5 


37.7 


14 


175 


35. 


15 


182 


36.4 


16 


203 


40. 6 


17 


-210 


42. 


18 


-182 


36.4 



Production well oil rates during fill-up equal 20 STB/Day. 
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FAULT BLOCK 4 
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FAULT BLOCK 5 




O Production Well 

Figure 13 Well Locations for Fault Blocks 3-5. 



TABLE VII 



WELL DATA— FAULT BLOCK 4 
Refer to Figure 13 



Well 

No. 


Rate, STB/Day 
(+Inj, -Prod) 


Thickness 

Ft. 


1 


88 


44 


2 


-100 


50 


3 


90 


45 


4 


-100 


50 


5 


86 


43 


6 


-94 


47 


7 


48 


48 


8 


-48 


48 


9 


-48 


48 


10 


48 


48 



Production well oil rates during fill-up equal 10% of initial rate. 



TABLE VIII 



WELL DATA— FAULT BLOCK 5 
Refer to Figure 13 



Well 


Rate, STB/ Day 


Thickness 


No. 


(+Inj, -Prod) 


Ft. 


1 


-50 


50 


2 


* 


50 


3 


-50 


50 


4 


* 


48 


5 


-50 


50 


6 




50 


7 


-46 


46 


8 


* 


46 


Production well oil rates during fill-up equal 3 STB/Day 



* Program calculates injection well rates based on reservoir thicknes 
and steady state material balance. 
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TABLE IX 

WELL DATA— FAULT BLOCK 6 



Refer to Figure 14 



Well 

No. 


Rate, STB/Day 
(+Inj, -Prod) 


Thickness 

Ft. 


1 


260 


56 


2 


240 


52 


3 


245 


53 


4 


245 


53 


5 


225 


49 


6 


-340 


58 


7 


-300 


60 


8 


-300 


60 


9 


-300 


60 


10 


-285 


57 


11 


-255 


51 


12 


-275 


55 


13 


300 


60 


14 


-300 


60 


15 


-270 


54 


16 


250 


48 


17 


305 


52 


18 


-275 


55 



Production well oil rates during fill-up equal 25 STB/Day. 
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4 



FAULT BLOCK 7 





Figure 14 Well Locations for Fault Blocks 6-8. 
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TABLE X 

WELL DATA— FAULT BLOCK 7 
Refer to Figure 14 



Well 

No. 


Rate , STB /Day 
(+Inj, -Prod) 


Thickness 

Ft. 


Well 

No. 


Rate , STB /Day 
(+Inj , -Prod) 


Thickness 

Ft. 


1 


234 


39 


13 


-282 


47 


2 


-264 


44 


14 


258 


43 


3 


234 


39 


15 


-288 


48 


4 


-258 


43 


16 


276 


46 


5 


234 


39 


17 


-288 


48 


6 


-258 


43 


18 


270 


45 


7 


228 


38 


19 


-282 


47 


8 


-240 


40 


20 


264 


44 


9 


222 


37 


21 


-276 


46 


10 


-228 


38 


22 


258 


43 


11 


-270 


45 


23 


-270 


45 


12 


252 


42 


24 


258 


43 



Production well rates during fill-up equal 25 STB/Day 



TABLE XI 



WELL DATA— FAULT BLOCK 8 



Refer to Figure 14 



Well 

No. 


Rate, STB /Day 
(+Inj, -Prod) 


Thicki 

Ft, 


1 


124 


31 


2 


-124 


31 


3 


124 


30 


4 


-124 


31 


5 


-168 


42 


6 


168 


42 


7 


-168 


42 


8 


160 


40 


9 


-168 


42 


10 


156 


39 



Production well oil rates during fill-up equal 15 STB/Day. 
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TABLE XII 

WELL DAT A — F AU LT BLOCK 9 



Refer to Figure 15 



Well 

No. 


Rate, STB /Day 
(+Inj, -Prod) 


Thicknes s 
Ft. 


Well 

No. 


Rate, STB/Day 
(+Inj , -Prod) 


Thicknes s 
Ft. 


1 


A 


27. 7 


14 


-150 


40 


2 


3jC 


30.4 


15 




40 


3 


-120 


29. 


16 


-150 


40 


4 


-160 


31.7 


17 


-200 


40 


5 


'I' 


30.4 


18 




40 


6 




35. 6 


19 


-150 


40 


7 


-135 


33. 


20 


A 


40 


8 


-200 


40. 


21 


-150 


40 


9 


❖ 


35.6 


22 


-150 


40 


10 




40 


23 




40 


11 


-155 


37. 6 


24 


V*/ 

V 


40 


12 


-200 


40. 


25 


-150 


40 


13 


•t* 

•T* 


40 









Production well oil rates during fill-up equal 15 STB/Day. 

❖ Program calculates injection well rates based on reservoir thickness 
and steady state material balance. 



\ 
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FAULT BLOCK 9 





X Injection Well 
O Production Well 



Figure 15 Well Locations for Fault Blocks 9 and 10. 
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TABLE XIII 

WELL DATA— FAULT BLOCK 10 
Refer to Figure 15 



Well 

No. 


Rate, STB/Day 
(+Inj, -Prod) 


Thickness 

Ft. 


Well 

No. 


Rate , STB/Day 
(+Inj, -Prod) 


Thicknes s 
Ft. 


1 


-76 


38 


21 




40 


2 




35 


22 


-80 


40 


3 


•f* 


40 


23 


-80 


40 


4 


-76 


38 


24 


-80 


40 


5 


-70 


35 


25 


* 


40 


6 


-80 


40 


26 


-76 


38 


7 


-74 


37 


27 




35 


8 




35 


28 




40 


9 


❖ 


40 


29 


-80 


40 


10 


-80 


40 


30 


-78 


39 


11 


-74 


37 


31 


-76 


38 


12 




34 


32 


-72 


36 


13 




40 


33 


-70 


35 


14 


-150 


40 


34 


-t" 


38 


15 


-150 


40 


35 


-76 


38 


16 




40 


36 


* 


37 


17 


-80 


40 


37 




36 


18 


-76 


38 


38 


* 


33 


19 


-72 


36 


39 


-66 


33 


20 


'i" 


34 









Production well oil rates during fill-up equal 8 STB/Day. 

❖ Program calculates injection well rates based on reservoir thickness 
and steady state material balance. 



APPENDIX B 
RATE ADJUSTMENT 



51 



ADJUSTMENTS OF INJECTION AND PRODUCTION RATES TO 
ACHIEVE UNIFORM WATER BREAKTHROUGH 
IN ALL FAULT BLOCKS 

When the model response was inspected section by section it was 
found that some of the sections were watering-out faster than others. 

This situation is also possible in the field. If the injection and production 
rates cannot be sufficiently adjusted, certain sections would be watered- 
out and shut-in while others would be run longer until ultimate recovery. 

If the rates can be adjusted enough the field can be operated so 
that all of the sections water -out at about the same time. The well rates 
used in the model were adjusted, with the total field injection rate held 
constant, to determine the effect of uniform water breakthrough from 
section to section. 

The results, as shown in Figure 16, indicate that for the same 
period of operation the oil production was increased and the water 



production was decreased. 



40h 
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Figure 16 Oil and Water Production With Well Rates Adjusted for Uniform Water 
Breakthrough Times. Predictions Made Using Streamline Model. 



APPENDIX C 



COMPUTER PROGRAM 

1. Program Listing 

2. Sample Data Deck 



*0 00 00 00*0 00 0*00 00 00 0000000 00 OOO OO OOOOOOOO O 0*0 0*0 oooooooo oooooo oo 
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PROGRAM OLBNKVTdNPUT, OUTPUT) ft l 

ft 2 
ft 3 

A STREAMLINE MODEL FOR WATERFLOOD SIMULATION A 4 

ft S' 
A 6' 

THIS MODEL IS SPECIFICALLY DESIGNED TO SIMULATE A ft 7 

WATERFLOOD OPERATION IN WhICH An OIL BAN* IS FORMED AND A 8 

MOVES THROUGH THE: RESERVOIR. ThE MODEL IS CAPABLE OF A 9 

HANDLING RESERVOIRS WITH THICKNESS VARIATIONS ftNO A 10 

IRREGULAR BOUNDARIES. ft 11 

4 12 

REFERENCES! A 13 

A 14 

LEBLANC. J.L. » AND CAUDLE, 8. H., A STREAMLINE MODEL FOR A 15' 

SECONDARY RECOVERY* SOC. PeT. EngR. J. (1971). 11*7. A is. 

RUST, C. 8,. A STREAMLINE MODEL FOR OIL BANK 0UILOJP IN A 17 

a waterflood, m.s. thesis, the university of texas at a 19 

AUSTIN (MAY. 1972). ‘ A 19 

WESSELS, J.W. , APPLICATION OF THE STREAMLINE RESERVOIR A 2fl 

MODEL TO RESERVOIRS HAVING VARIATIONS IN FORMATION A 2l 

THICKNESS. M.S. THESIS. THE UNlYERSIlY OF TEXAS AT AJSTIN A 22 

(MAY, 1973). A 23 

A 24 

DOCUMENTED BY! K.E.GOLTZ (MARCH, 1975) A 25' 

A 2S 
A 27 

»•••«•••««•<>«••«••••»•*»»»»•••«•»• A 23 

* SCHEMATIC PROGRAM flow CHART * A 23 

a 3J 

A 31 
A 32 

... A 33 

I read input data i a 34 . 

A 35 

I A 3 S' 

A 37 

I COMPUTE STREAMLINE STARTIN3 POINTS. I A 33 

I SOURCE AND SINK PRESSURES I A 39 

.......................................... A 40 

I A 41 

................................ A 4 2 

1 ITERATE Each TIME PERIOD I A 43 

............................... A 44 

I A 4 5' 

................................... a 4 S' 

— ..........I ITERATE each PRODUCTION WELL I A 47 

I A 43 

I I A 49 

I ........... ................ A 50 

I -j iterate EACH STREAMLINE I A 51 

I I - — A 52 

II I A 53 

I I .................. a 54 

I I ——I COMPUTE VELOCITIES AND TIME INCREMENTS I A 55' 

III I AT EACH POINT ALONG STREAMLINE I A 5S 

I I I A 57 

III I A 59 

I I I — A 59 

III I MOVE FLUID PARTICLE ALONG STREAMLINE USING I A SO 

III I FINITE DIFFERENCE ApPROX. OF DARCTS LAW I A 61 



0*0*00 0*00 O O O O O OO 0*0 O O 0*0 OO OO OOOOOOO O O O O O O O O 0*0*0 O O O O O 0*0 0*0*0 0*00 '0 0*0 o o oo 
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I I I — — — ft 62 

III I ft 63 

I I I — * — — — ft 6ft 

III I ACCUMULATE STREAMLINE PRODUCTION I A 65' 

I I I — -- A 65 

III I A 67 

I I I A 69 

III I TEST FOR BREAKTHROUGH OF FLUID I ft 69 

III I PARTICLES INTO PRODUCTION WELLS, I ft 70 

III I ADJUST STREAMLINE PRODUCTION I ft 71 

III I RATES AS REQUIRED I A 72. 

I I I - — A 73 

III I A 7ft 

III NO — — A 73' 

I i ..... ..i TEST FOR END OF TIME PERIOD I A 76 

I I — ... ft 77 

II I YES ft 79 

I I NO ....... & 79 

i ....... ———.i test for all streamlines i a so 

I ... A 91 

I I YES A 92 

I NO ........ ... A 03 

.... — -..-..-.-I TEST for ALL PRODUCTION WELLS I A 9ft 

.... ......... ...... ft S3' 

I YES ft 95 

..... ....... ......... ft 97 

I PRINT/PLOT WElE and FIELD RESULITS I A 99 

......... A 99 

I ft «0 

NO — .... A 9i 

.....................I TEST FOR TOTAL time I A 92 

... A 93 

I YES A 9ft 

... A 95' 

I STOP I ft 96 

.......... A 97 

A 99- 
A 99 
A 100 

A 101 

INPUT DATA DEFINITIONS A 102 

A 1 0 3 

POR a POROSITY (DECIMAL) A loft 

SC« = SATURATION CONNATE water (DECIMAL) A 1 0 5' 

SOI * SATURATION OIL INITIAL (DECIMAL) A 106' 

SOR a SATURATION OIL RESIDUAL (DECIMAL) A 1 0 7 

SGR = SATURATION GAS RESIDUAL (DECIMAL) A 109 

sgor s solution gas-oil ratio (scf/bbl at res. conditions) a 109 

np c number of production wells a lio 

NI S number of injection WELLS a 111 

NBW e NUMBER OF IMAGE WELLS USED TO BOUND PATTERN ft 1 1 2 

NBP b NUMBER OF POINTS USED TO DEFINE BOUNDARY A 113 

NST * NUMBER OF STREAMLINES FOR LOWEST RATE PRODUCTION WELL ft lift. 

NXH e NUMBER OF X-DIRECTlON ENTRIES IN THICKNESS GRID ft 115* 

NYH a NUMBER OF Y-DIRECTION ENTRIES IN THICKNESS GRID ft 1 1 6 

NEOIT a o OR 1 ft 117 

IF 6 OR BLANK-PROGRAM WILL BALANCE PROD/INJ RATES A 119 

IF l PROGRAM INTERNAL BALANCE FEATURE IS DEFEATED A 119 

H ■ AVERAGE RESRVOIR THICKNESS, FT (USED FOR IMAGE WELLS) A l2o 

RI » RADIUS OF STARTING AND ENDING CIRClES.fT A 121 

DELT b TIME INCREMENT (DAYS) FOR UPDATING CONSTANT ORESSJRE A 122. 

CORRECTION FACTOR, SINCE DELT MUST START OUT SMALL BUT A 123 
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c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

t 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 



BECOMES LESS CRITICAL AS ThE SOLUTION PROGRESSES, 

the value o ? i oelt is wultipIed by te.\i after the 
OIL Bank BREAKS THROUGH in The FIRST STREAMLINE. 

OELP ■ TIME INCREMENT FOR PRINT AnD PLOT (DAYS) 

FIRST PRINT/PLOT IS MADE AT CELT, THEN EVERY DELP. 
TMX s MAXIMUM'TIME FOR SOLUTION (DAYS) 

QMNP = LOWEST PROOJCTION WELL RATE (STB/DAY) 

OMNI s LOWEST INJECTION WELL RATE (ST3/DAY/FT) 

xmx * absolute x dimension at which streamlines arei cut off- 

YMX o ABSOLUTE Y DIMENSION at WHICH STREAMLINES arei CUT OFFi 

GM8.OM3.WMB » GAS. OIL AND WftTER MOBILITIES 

XW ( I ) . YW ( 1 ) b X.Y COORDINATES OF WELL LOCATION 

HI(I) s. THICKNESS OF RESERVOIR AT WELL. FT 

Q(I> * WELL' INITIAL' OIL PROOJCTION RATE ( ST 3/OAY ) 

QO(I) * WELL PRESENT (FILL-UP) OIL PRODUCTION RATE (STB/OAY) 
GOR(I) * WELL GAS-OIL RATIO (SCF/SBL AT RES. CONDITIONS) 

X8 ( I ) . Y8 ( I ) s X.Y COORDINATES OF BOUNDARY POINTS 

xh(I) ,YH( i) = x.y Coordinate values for thickness grid 

HTd.J) = VALUES OFi THICKNESS CORRESPONDING TO COORDINATES 
OF THICKNESS GRIO 






COMMON/AAA/ XW(125) .YW(125) , Hi ( 1 25) .0 ( 1 25) . XB < 125) . 

1 YB ( 1 25) .NBP.NP.NI.NBW.H 
COMMON/BBB/ C (55.55) ,D( 55) .01 (5S) 

COMMON/CCC/ X (80*20) ,Y<80»20> *X0 (80.20) . Yu (90.20 > . 

1 WQ(80. 20) .SCL. XMX. YMX.NST.OMNP.RGO <BO»20 ) * 
2SG0R.NRIW.TMEI 

DIMENSION TSUM(80.20) .PP<80.20> .Pi <80*20> »C0<80.20) * 
ICW(B0.20) .00(125) .50R( 125) .00(80.20) . 

2XH(20> . YH ( 20 ) .HT (20,20) 

READ 530, POR.SCW.SOI.SOR.SGR.SGOR 
np.ni.nbw.nbp.nst 
nxh.nyh.nedit 

H.R I, OELT, DELP. TMX, QMNP 
OMNI »XMX, YMX,GM8»0MB,wMB 



READ 

READ 

READ 

READ 



540, 

540, 

530, 

530, 

DELTTbDELT*10. 

NC0UNT«0 

PLTFTb0.9«XMX/YMX*0,45 

TOTFTso.7*PLTFT 

NPlW=NP*Nl 

NPIBwbNPIw*NBW 

READ 550, <XW(I) ,YW(I) ,HI (I) ,0(1) ,QO(I) ,GOR(I) ,I=1,NPIBW) 
READ 560, ( X 8 ( I > ,Y8(I) ,1 = 1. N3P) 

CONVERT SURFACE RATES TO RESERVOIR RATES 
FVFi b 1.2 



DO 9 1 ■ 1 , NPI W 
Q(I)»Q(I)*1»2 
QO(I)bOO(I)* 1.2 
QMNPbqmNPM. 2 
QMNl*OMNI*l . 2 
PRINT 590 
PRINT 480. 

PRINT 49o. 

PRINT 500. 

PRINT 510. 

PRINT 520. 

PRINT 550. 

PRINT 530. 



NP.NI.N3W.NBP.NST 
H.DELT. TMX.Rl 
DELP. XMX, YMX.QMNP 
POR.SCW.SOI.SOR.SGR 
OMB.WMB.GMB 

(XW(I> , YW ( I ) , HI (I) . Q ( I ) ,00(1) . GOR ( I ) »I*1,NPIBW) 
( XB ( I ) , YB ( I ) ,1*1, NBP) 



A 


m 


A 


125 ' 


A 


126 - 


A 


127 


A 


128 - 


A 


129 


A 


130 


A 


131 


A 


132 


A 


133 


A 


13 V 


A 


135 ' 


A 


13 S 


A 


137 


A 


138 


A 


139 


A 


140 


A 


Hi 


A 


H? 


A 


143 


A 


144 


A 


H 5 


A 


14 S 


A 


147 


A 


149 


A 


149 


A 


150 


A 


151 


A 


152 


A 


153 


A 


154 


A 


IBS' 


A 


1 5 S' 


A 


157 


A 


159 


A 


159 


A 


160 


A 


161 


A 


162 


* A *1 6 3 


* A * 1 64 


A 


165 ' 


A 


166 


A 


167 


A 


169 


A 


169 


A 


170 


A 


171 


A 


172 


A 


173 


A 


174 


A 


175 


A 


176 


A 


177 


A 


179 


A 


179 


A 


1 9 0 


A 


1 B 1 


A 


152 


A 


1 B 3 


A 


184 


A 


IBS’ 



57 



C A 13V 

C BALANCE INJECTION WELL RATES CONSIDERING RESERVOIR A 187 

C THICKNESS AND STEADY STATE OPERATION A 193. 

t A 199 

IF(NEDIT ,GT. 0) GO' TO 31 A 1 90 

SMQUo.O ' A 191 

SMQP*o«0 A 192 

DO 20 1*1 * NPIW A 193 

IF (Q(I).LT.O.O) GO TO 10 A 194 

SMQI»SMQI*Q(I) A 195' 

GO TO 20 A 195. 

10 CONTINUE A 197 

SMQPeSMQP*Q(I)/HI (I) A 199 

20 CONTINUE A 199 

DO 30 1*1 * NPIW A 200 

IF (0(1) .LT. 0.0) GO TO 3o A 201 

Q(I)=-a(I)*SMQP*HI (D/SMQI A 202 

30 CONTINUE A 203 

GO TO 32 A 204- 

31 PRINT 595 A 203' 

32 CONTINUE A 205 

C A 20 7 

C check if oil bank will forh a 209 

C A 20 9 

IF (SOR.LT.SOI) GO TO 40 A 21 0 

PRINT 570 * A 211 

GO TO v70 A 212 

♦0 CONTINUE A 213 

C A 2U 

C CONVERT field UNITS TO CGS A 215 

C A 215 

DO So I=1.NPI6W A 217 

XW(I)bXW(I)*30.48 A 219 

YW(I)»YW(I)«*3o.49 A 219 

00(I)=00(I)*1.8401 A 220 

HI ( I ) *Hl < I ) *30 »*8 A 221 

Q(I>*Q(I)*1.840l A 22? 

IF ( I .LE.NPIW) 30 TO So A 223 

Hill) bH*30»48 A 224 

50 CONTINUE " A 225' 

DO 60 1*1 • N8P A 225 

XB(I)bXB(I)*30.49 A 227 

YB(I)bYB(I)*30.A8 A 223 

60 CONTINUE A 229 

HbH*3O.A0 A 230 

DO 70 1*1. NXH A 231 

READ 580. (HT(I.J) »J*1.NYH) A 232 

70 CONTINUE A 233 

READ 530. (XH(I) .1*1. NXH) A 234- 

READ 580. ( YH ( I ) .1*1 . NYH) A 235' 

DO 80 I*1.NXH A 235- 

DO Go J*1 . NYH A 237 

HT(I»J)*HT (I.J)*30.48 A 239 

80 CONTINUE * A 239 

DO 90 1*1 . NXH A 240 

XH(I)bXH(I)*30.48 A 241 

90 CONTINUE A 242 

DO 100 1*1. NYH A 243 

YH(I)»YH(I)*30.48 A 244 

100 CONTINUE A 245' 

R1*RI*30.48 A 246 

DELT«DELT*86400. A 247 



O'O OO O 00*0*000 O O O O 00«00'00>0 OO O OO O OO O O OO OO 00*0*0 O O OO'OO 
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OEUTT«DELTT*B6400. A 249 

DELP»DELP*86400, A 249 

TMXsTMX«*e6400. A 25o 

QMNP«OMNP*1.840l A 251 

XMX«XMX*30.4B a 252. 

YMXaYMX*30 .48 A 253 

QMNI=OMNI*l.B40l A 254 

A 255' 

SUBROUTINE bound CALCULATES The NECESARY RATES A 255. 

FOR PREVIOUSLY SELECTED IMAGE WeLLS SO THAT THE A 257 

STREAMLINES ARE CONFINED WITHIN THE RESERVOIR ‘ A 259 

BOUNDARIES. ' A 259' 

A 260 

CALL BOUND A 261 

PRINT 600 A 262 

A 263 

CHECK OMNI AND CALCULATE VELOCITIES NEAR' A 264 

INJECTION AND PRODUCTION WELLS A 265' 

A 266' 

00 1 1 6 I *1 1 NPI 8W A 267 

IF (Q(I).LT.O.O) GO TO 110 A 269 

IF (Q(I)/MI (I) .LT.3MNI) QMN I «Q ( I ) /HI ( I > A 269 

10 CONTINUE A 270 

VLCI*o.5*»QMNI/RI A 271 

VLCPr 0 ,4*QMNP/H/RI A 272 

NTSL*0 A 273 

A 274. 
A 27 5' 

THIS SECTION STARTS A STREAMLINE ON THE CIRCLE OF 4 A 275 

PRODUCTION WELL AND FOLLOWS IT BACKWARDS TO LOCATE THE A 2.77 

PROPER INJECTION WELL. This PROCESS is CARRIED OUT for a 279 

ALL PRODUCTION WELLS AND ALL STREAMLINES. THE A 273 

STREAMLINE ST0P5 AT A POINT WIThIN THE CIRCLE OF AN A 290 

INJECTION WELL. THIS POINT BECOMES THE STARTING POINT A 291 

WHEN The STREAMLINE IS FOLLOWED in the FORWARO OIRECTION a 282, 

TO ACCUMULATE THE. PRODUCTION HISTORY. A 293 

AS THE STREAMLINES ARE BEING FOLLOWED A PLOT OFi THEIR *A*284 

TRACE is made. (NOTE THaT ALL PLOT STATEMENTS HAVE *A*283' 

ASTERISKS In THEIR STATEMENT NUMBERS.) *4*296' 

A 287 
A 289 
A 283 

PROGRAM VARIABLES A 290 

A 291 
A 292. 

J » INDEX FOR WELLS A 293 

k « index for streamlines a 294 . 

A 295' 

XIiYI « POSITION ON STREAMLINE OF FLUID PARTICLE BEUNG TRACED A 295. 

VX » VY B VELOCITY OFi FLUID PARTICLE IN X AnD Y DIRECTIONS A 297 

PP « PRESSURE AT PRODUCTION WELL (SINK) A 299 

PI * PRESSURE AT INJECTION WELL (SOURCE) A 299 

X(JtK) »Y(J.K) a STARTING POINTS Fo« WATER FRONT FROM INJ. WELL A 3oo 

XO (J»K) fYO(J.K) a STARTING POINTS FOR OIL FRONT FROM INJ^ WELL A 301 

TSUM(J.K) a SUMMATION OF TIME A 302 

RGO(J.K) » STREAMLINE GAS-OIL RATIO A 303 

OQ(J.K) a STREAMLINE OIL RATE A 304 

CO(J.K) a STREAMLINE CUMMULATIVE OIL PRODUCTION A 305' 

WQ(J.K) a STREAMLINE WATER RATE A 306. 

CW(J.K) s STREAMLINE CUMMULATIVE WATER PRODUCTION A 307 

SCL b PLOT SCALE FACTOR *A**3o9 

NSL « NUMBER OF STREAMLINES FOR ThE PRODUCTION WELII A 3o3 
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c 


BEING WORKED ON 




A 310 


c 


DT * time increment for small step taken, step size 


IS 


A 311 


c 


EQUAL TO .25**1. RI NEED NOT BE SMALLER THAN 




A 312, 


c 


1/50TH OF THE DISTANCE BETWEEN WELLS. 




A 313 


c 






A 31 A 








A 315 


c 






A 315 


c 






A 317 




SCL S 5. 238/ YMX 




*A»3l 9 




XM*2.25*XMX*SCL 




* A*3l 9 




CALL 63NPLT (ALPLOT.XM.25. 25) 




*A*32o 




CALL PLT (1.125*XMX*SCL*l.l25*YMX*SCL*"3) 




*A*32i 




CALL PLT (XS(NBP)*SCL*YB(NBp)*SCL,3) 




*A*322 




DO 120 1*1 * N3P 




•A* 323 




CALL PLT <X8(I)*SCL4YB(I)*SCL*2) 




*A*32a 


120 


CONTINUE 




*A*325' 




DO Uo I *1 * NPI W 




*A*325, 




IF (Q(I) .SE.O.O) GO TO 130 




«A*3?7 




CALL SYMBOL (XW ( I ) »SCL * Y W ( I ) *SCL * . 07 * 1 ♦ 0 »-l > 




*A*323 




GO TO 1 AO 




*A*329 


130 


CONTINUE 




*A*33o 




CALL SYMBOL (XW(I)»SCL*YW(I)*SCL*.07*T.0t-l) 




*A*33l 


140 


continue 




*A*332 




DO 220 Usl.NPIW 




A 333 




IF (Q(J) .GE.O.O) GO TO 220 




A 334 




NSL=-NST*Q(U)/3MNP*0.5 




A 335 


C 






A 335 


c 


THE STREAMLINE STARTING POSITIONS ARE EVENLY SPACED 


A 337 


c 


AROUND THE PRODCTION WELL 




A 339 


c 






A 339 




DO 21 o Ksl.NSL 




A 340 




XI*XW(U) *RI*COS(l .0*K«6.2832/NSL) 




A 341 




Y1«YW( J)*RI*SIN(1.0*K*6.2B32/NSL) 




A 34? 




CALL PLT (XI*SCL.YI*SCL.3) 




* A*343 




PP( J.K)=0.0 




A 344 




DO 150 L*1.NPI8W 




A 345' 




PP( J.K)=PP(J,K)-Q(U)/HI (L)*ALOG( (Xl-XW(L) )**2*(YI-YW(L> )**2) 


A 345 


150 


CONTINUE 




A 347 


160 


continue 




A 349 




IF (ABS (XI) .GT.XMX) GO To 180 




A 349 




IF (ABS(YI) .GT.YMX) GO To 180 




A 350 




VX«0.0 




A 351 




VY*0 . 0 




A 352 




DO 17o L= 1 . NP I BW 




A 353 




DEN*(XI-XW(L) )**2* (YI-YW(L) ) **2 




A 354 




VXsVX-Q(L)*(XI-XW(L) )/(DEN*Hl <L) ) 




A 355' 




VY*VY-3(L)*(YI-YW(L) )/(DEN*Hl <L) ) 




A 355 


170 


continue 




A 357 




VT*SQRT (VX**2*VY**2.) 




A 359 




DT«0.25*RI/VT 




A 359 




xi*xiot*vx 




A 360 




YI«YIOT*VY 




A 351 




CALL PLT (XI*SCL.YI*SCL*2) 




*A*352 


C 






A 363 


c 


CHECK IF THE STREAMLINE HAS REACHED AN INJECTION 


WELL. 


A 364 


c 


the velocity of the fluid particle is checked first to 


A 365 


c 


SEE If IT IS LARGE ENOUGH to be in The VICINITY OF 


AN 


A 365 


c 


INJECTION WELL. IF IT IS, EACH INJECTION WELL IS 


then checked. 


A 357 


c 






A 369. 




IF (VT.LE.VLCI) GO TO 16o 




A 359 




DO 200 LL*1«NPIBW 




A 370 




IF (Q(LL) .LE.0.0) GO TO 200 




A 371 



60 



RAD«SQRT( (XI-XW(LL) ) **2» ( Vl-YW (LUj )**2> A 372 

IF (RAO.GT.RI) 30 TO 200 A 373 

180 CONTINUE A 371 

NTSL=NTSL*1 A 375' 

X ( J» K ) cX I A 37S 

Y(J*K)*YI A 377 

TSUM(J.K) cOtO A 373 

XOU.KHXI ' A 37? 

YO ( J» K) °YI A 38o 

PI(J.K)»0.0 A 381 

00 190 LK=1»NPIBW A 382 

PI (J.KHPI (J.K)-Q(LK)/HI (LM*ALOG( (X I-XW (LK> ) **2* ( YI -YW (LK) >**2) A 38 3 

190 CONTINUE A 33 A’ 

RGO ( J, K) =GOR < J) A 383 1 

IF (SCW*S0I*SGR.GE.1.0) RGO ( J.K) »S<SOR A 38$ 

OQ(JiK)=QO(J)/NSU A 397 

CO(J*K)=0.0 A 399 

WQ(J,K)«0.0 A 39? 

GO TO 210 A 39o 

200 CONTINUE A 391 

GO TO 160 A 392 

210 CONTINUE A 393 

CW(J.K)=0.0 A 39V 

220 CONTINUE A 395 

CALL PUT (O.OtO.O* 999) «A*3?s 

PRINT 610 ' A 397 

C A 393 

C END OF RUNBACK SECTION A 399 

C A AOD 

CTR»1.0 A A01 

PRINT 590 A »0? 

TSTOPsDELT A A03 

TSMP=OELT A AOA 

NBP1=N3P*1 A A05 1 

230 CONTINUE A Aos 

SMCTYeO.O A AO 7 

C A *09 

C A Ao? 

C THIS SECTION STARTS ThE STREAMLINE AT THE INJECTION A A-lo 

c well found above ano follows it in the forward direction a ah 

C TOWARD THE PRODUCTION well. the MOVEMENTS of THE' OIL AND A Al? 

C WATER fronts ARECLOSELY WATCHED TO DETERMINE WHEN THEY A An 

C BREAKTHROUGH INTO THE PRODUCTION WELL. THE PRODUCTION A A1a 

C HISTORIES ARE ACCUMULATED AS THE FRONTS MOVE. A A 1 5' 

C LOOPS A20 AND HO ARE DONE FOR ALL WELLS AND ALL A Al$. 

C STREAMLINES UNTIL TIME EQUALS TsTOP. THEN THE CONSTANT A V 1 7 

C PRESSURE CORRECTION FACTOR <CTR) IS UPDATED AND THE A An 

C RESULTS ARE PRINTED IF DESIRED. a H? 

C A A 2 j 

c A A 21 

C A A?2 

C ADDITIONAL' PROGRAM VARIABLES A *23 

C A A2v 

C CTY E CONDUCTIVITY RATIO A A25' 

C CTR » CONDUCTIVITY RATIO CORRECTION FACTOR A *2$ 

C TSTOP « TIME TO STOP AND UPDATE CTR A »27 

C TSMP e TIME TO PRINT ANO PLOT A *29 

c xitYi e location of. water front on streamline being followed a n? 

C X01.Y0I « LOCATION OF OIL FRONT ON STREAMLINE BEING FOLLOWED A *30 

C A A3i 

C ***************************************************************** A A 3 2 

C A A33 
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DO 426 J*1 iNPIW 


A 


*3*. 




IF ( Q ( J) .6E.0.0) GO TO 420 


A 


435' 




N$L*-NST»Q( 0)/QMnP*O,5 


A 


436 




DO 410 K=1,NSL 


A 


437 




CTY»0. 0 


A 


439- 




IF (ABS (X < J,K) ) .GT.XMX) GO TO 4l0 


A 


439 




IF (ABS(Y(J,K) ).GT,YMX) GO TO 4l0 


A 


449 




IF (WQ(J,K> .GT.0.0) GO TO 400 


A 


441 




X I»X K) 


A 


442. 




YI=Y(J,K) 


A 


443 




TSUMI»TSUM( J f K> 


A 


444 




XOI*XO(J*K) 


A 


445' 




YOI«YO(J*K) 


A 


446 


240 


CONTINUE 


A 


447 


c 




A 


449 


c 


determine thickness OF reservoir at location of FRONTS 


A 


449 


c 




A 


450 




DO 250 MsltNXH 


A 


451 




IF (XI.UT.XH(M> ) GO TO 260 


A 


452 


250 


CONTINUE 


A 


453 




MaNXH 


A 


454 


260 


CONTINUE 


A 


455’ 




DO 27o Nal » NYH 


A 


456 




IF (YI.LT.YH(N) ) GO TO 280 


A 


457 


270 


CONTINUE 


A 


459 




NaNYH 


A 


459 


280 


CONTINUE 


A 


460 




htp=ht(M,n» 


A 


461 




IF (RGO(UfK) .LE.SGOR) GO TO 360 


A 


462 




DO 29o Msl.NXH 


A 


463 




IF (XOI.LT.XH(M) ) GO TO 300 


A 


464- 


290 


continue 


A 


465' 




MsNXH 


A 


466 


300 


continue 


A 


467 




00 310 N«1,NYH 


A 


469 




IF (YOI.LT.YH(N) ) GO TO 320 


A 


469 


310 


continue 


A 


47 0 




NaNYH 


A 


471 


320 


CONTINUE 


A 


472 




HTPO«HT(M»N) 


A 


473 


C 




A 


47 4 


c 


THIS SECTION MOVES THE OIL AND WATER FRONTS AND ACCUMULATES 


A 


475 1 


e 


PRODUCTION BEFORE. THE OIL FRONT REACHES A PRODUCTION WELL, 


A 


476 


c 




A 


477 




vxo*o,o 


A 


479 




VY0=0,0 


A 


479 




P0°0.0 


A 


49o 




VXbO.O 


A 


491 




VY«0,0 


A 


492 




PaO. 0 


A 


493 




DO 336 Lal'NRIBW 


A 


494- 




SQO= (XOI-XW (L) > **2* (YOI-YW (L> >**2 


A 


405 




P0«P0-3(L>/HI (L) *ALOG(SQO) 


A 


406 




VXOaVXO*Q(L)/HI (L) » (XOI-XW (L> )/SQO 


A 


497 




VYO«VY0*Q(L> /HI <L> » (YOI-YW (L> )/SQO 


A 


499 




SQ» (XI -Xw ( L) )*»2* (YI-YW(L) )»*2 


A 


499 




P»P-Q(L) /HI (L) *AL03(SQ) 


A 


490 




VXaVX»3(L)/HI (L> • ( XI-XW (L) ) /SO 


A 


491 




VY«VY*a(L)/HI(L)*(YI-YW(L) )/SQ 


A 


492 


330 


continue 


A 


49 3 




IF (PO.LT.PP(U.K) ) PO»PP(J*K) 


A 


494- 




CTY«(PP( J»K) -PI ( J*KJ )/( ( (PP( J*K) -POI/0MB* (PO-P>/OM3* (P-PI ( J,K> > / 


A 


495' 



0 0*00 O OO O 0 0<0<0<0 0 0 
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1 WMB)*GMB) a *96 

VIX=VX/(6.2832«P0R | *(1-SCW-S09-SGR) ) #CTY/CTR*Hl < J) /HTP A *97 

VLY»VY/(6.2832*POR*(l-SCW-SoR-SGR> ) *CTY/CTR*Hl ( J) /HTP A *99. 

V£L=SQRT ( VLX**2*VLY**2) A *99 

DT«RI/VEl A 500 

CRF=l,o*OQ( J.K)*NSL*CTR/(Q< J) *CTY> A 901 

VLX0sVX0/(6.2832*P0R*(1-SC*'-S0I-SGR) ) *CTY/CTR*CRF*HI ( JJ/HTPO A 902 

VLY0=VY0/(6;2832*PDR*tl-SC«»-S0l-SGR) )*CTY/CTR<*CRF*HI (JJ/HTPO A 903 

VE10»SQRT(VUX0**2*YLY0»*2) " A 90* 

IF (VELO.GT.VED DT=RI/vELO A 503' 

IF (TSUMI*DT .GT.TSTOP) DT*TsTOP-TsUMI A 906 

XIbXI*DT*VLX ' A 907 

YI B YI*3T*VLY A 909 

TSOMI=TSUMI»OT A 909 

X0I=X0I*0T*VLX0 A 9 1 o 

YOI»YOI*OT*VLYO A 511 

A 91 ?. 

CHECK IF OIL FRONT HAS REACHEO A PRODUCTION WELL* A 513 

THE WATER FRONT IS ALSO CHECKED JUST IN CASE IT H4S GOTTEN A 5l * 

AHEAD OF'TrtE OIL FRONT. THIS OF COURSE IS NOT NORMAL A 515' 

BUT CAN HAPPEN I Fi THE OIL PRODUCTION RATE DURING A 5l 6 

FILL-UP IS TO LARGE. A 517 

A 9ls 

IF (SQRT((XI-XW(J))»*2*(YI-YW(J))*»2).LE,RI) SO TO 390 A 319 

CO(J.K)=CO< J.K)*DT»OQ(J.K) A 92j 

IF ( VELD .LE. VLCP) GO TO 35o A 521 

IF(SQRT((X0I-XW(J))*<*2*(YOl-YW(J))«*2).LE.RI) GO TO. 3*8 A 329 

DO 3*9 LX=1.NPI8W A 523 

IF ( U ( LX) .GE. 0.0) GO TO 3*9 A 52* 

IF(SQRT((X0I-XWfLX))**2*(Y0l-YW(LX))*«2).LE.RD GO TO 3*B A 525' 

3*9 CONTINUE A 526 

A 527 

WHEN THE OIL FRONT REACHES A PRODUCTION WELL THE A 529- 

STREAMLINE GAS-OIL RATIO IS SET EQUAL TO SGOR, A 529 

A 530 

00 TO 350 A 531 

3*8 RGO(J,K)=SGOR A 532 

NCOUNT ®1 A 333 

350 CONTINUE A S3* 

IF (TSUMI.LT. TSTOP) GO TO 2*0 A 535 

X(J»K)oXI A 336 

Y(J»K)«YI A 537 

XO(J»K)=XOI A 539 

YO(J»K)*YOI A 539 

TSUM( j,K)cTSUMI A 5 *o 

GO TO *10 A 5*1 

A 5*2 

THIS SECTION MOVES THE WATER FRONT AND ACCUMULATES A 5*3 

the production AFTER OIL FRONT BREAKTHROUGH A 5**. 

A 5*5 

360 CONTINUE A 5*5 

VX»0,0 A 5*7 

VYsO.0 A 5*9 

PeO.O A 5*9 

00 37o L«=1,NPI8W A 5S0 

SQ»(XI-XW(L) )**2*( YI-YW(L) )*»2 A 551 

P*P-G(L) /HI (L) ®AL03(SQ) A 352 

VX»VX*Q(L)/HI (L)*(XI-XW(L) )/SQ A 553 

VY»VYO(L)/HI<L)*<YI-YW(L))/SQ A 55* 

370 CONTINUE A 555 

IF (P.LT.PPI J.K) ) P«PP(J»KJ A 355' 

CTY«(PP( J.K)-PI ( J.X) )/( ( (PP( JfK)-p)/OMB. (P-RJ (J.K) ) /WMB) *GM8) A 557 
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VLX»VX/(6.2832#P0R» (1-SCW-S0R-SGR> >«CTY/CTR»Hl < J) /HTP 


A 


sss 




VLY c VY/(6.2832*P0R»(l-SCW-SOR-SGR) )*CTr/CT«»Hl (J)/HTP 


A 


559 




VELbSQRTIVLXbob^VLVbbE) 


A 


560 




ot»ri/vel 


A 


561 




IF (TSUMI*OT,GT,TSTOP) DT»TSTOP.TsUMI 


A 


562 




XIbXI*0T*VUX 


A 


^63 




YIoYI*DT*VLY 


A 


56*. 




TSUMIcTSUMI*OT 


A 


565’ 




OQ(J.K)=-Q(J)<*CTY/(CTR*NSL> 


A 


566 




CO(JtK)=CO(J.K)*OT*OQ( J.K) 


A 


567 


c 




A 


569 


c 


CHECK FOR WATER FRONT BREAKTHROUGH 


A 


569 


c 




A 


570 




IF (SQRT((XI-XW(J)J**2*(YI-YW<J))*°2).LE.RI) GO TO 390 


A 


571 


380 


CONTINUE 


A 


572, 




IF (TSUMI.LT.TSTOP) GO TO 2*0 


A 


573 




X ( Jt K ) bX I 


A 


57* 




Y(JtK) =YI 


A 


575 1 




TSUM( J,K)=TSUMI 


A 


576 




XO(J.K)=XOI 


A 


577 




YO( JtK) =YOI 


A 


579 




GO TO *10 


A 


579 


390 


continue 


A 


580 


c 




A 


581 


c 


THIS SECTION IS USED WHEN THE WATER FRONT FIRST 


A 


592 


c 


BREAKS THROUGH. 


A 


553 


c 




A 


58* 




RGO ( J. K) =SGOR 


A 


585 




CTY=WMB/GMB 


A 


586 




WQ(J.K)=-Q( J)«CTY/(NSL*CTR> 


A 


587 




OQ(J»K)=0.0 


A 


589 




X ( J » K ) a X I 


A 


599 




Y< JtK) sYI 


A 


590 




TSUM{J,K)=TSUMI 


A 


591 




XO ( J t K) s XOI 


A 


592 




YO( JtK)=YOI 


A 


593 




CW< JtK)»(TSTOP-TSU>1I)*WO(JtK) 


A 


59*. 




GO TO *10 


A 


593‘ 


400 


CONTINUE 


A 


596 


c 




A 


597 


c 


ACCUMULATE PRODUCTION AFTER WATER BREAKTHROUGH 


A 


599 


c 




A 


599 




CTYcWM3/GMB 


A 


500 




WO (JtK) =-Q(J)*CTY/(N5L*CTR) 


A 


501 




CW(JtK)=CW(JtK)*OELT*WQ(JtK) 


A 


602 


410 


smcty=$mcty*cty 


A 


603 


420 


continue 


A 


6o *• 




CTR»SMCTy/NTSL 


A 


605 




IF (TSTOP.LT.TSMP) GO TO *50 


A 


606 


C 




A 


607 


C 


AT SPECIFIED TIME INTERVALS 


A 


609 


c 


ALL STREAMLINES ARE SUMMED AND RESULTS ARE PRINTEO 


A 


609 


c 


FOR EACH WELL'ANO. FOR THE ENTIRE FIELD. 


A 


610 


c 


PLOTS ARE MADE Of i THE OIL FRONT AND WATER FRONT POSITIONS. 




c 




A 


612 




ORT*0 • 0 


A 


613 




WRT e 0 , 0 


A 


61*. 




SOIL»0.0 


A 


61 3‘ 




SWAT>0.0 


A 


616 




DO **0 JnitNPIW 


A 


617 




IF (Q(J).GE.O.O) GO TO **0 


A 


619 




GR°0. 0 


A 


619 
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0R»0.0 A S>2 0 

WR»0.0 A 621 

OPRD*0 • 0 A 622. 

WPRO»o.O A 623 

NSL 3 -NST*Q( J)/QMNP*0.5 A 62*. 

DO *30 Kel.NSL A 625' 

OPRD=OPRO*CO(J»K)/158990. A 625. 

WPRD»WPRD*CW (J*K)/15B99o* A 627 

GR 3 GR.RGO( j»K)*OQ< U.K) /l. 8*01 A 623 

OR»OR*OC)<J»K)/l .8*01 A 62? 

WR=WR*WQ(J*K)/\.8401 A 630 

*30 CONTINUE A 631 

GSR=GR/OR ft 632 

WOR=WR/OR ft 633 

EX«XW< J)/30.*8 ft 63* 

WY 3 YW ( J) /30 • *8 ft 635' 

PRINT 620 • U.EX » WY ft 635 

PRINT b3o» GR.OR.WR: ft 637 

PRINT 6*0 » OPRD.WPRD.GSR.WOR ft 633 

ORT=ORT *OR ft 63? 

WRT=WRT*WR ft 6*o 

SOILaSOlL*OPRD ft 6*1 

SWAT=SttAT**PRD ft 6*2 

**0 CONTINUE ft 6*3 

TMEI=TSTOP/86*00.0 A 6**. 

PRINT 650* ORT.WRT ft 6*5' 

PRINT 65l* SOIL*SWftT A 6*5 

PRINT 66o* TmEI.CTR A 647 

C *ft«6*3 

C MAXIMUM PLOT LENGTH IS 11 FT. *4*6*? 

C PLOT LENGTH IS ESTIMATED AND IF IT IS GREATER THAN *4*650 

C 11 FT* THE PLOTTING IS STOPPED. *ft*65l 

j *4*552 

IF ( ( 1 1 .0-T0TFT) .LT. PLTFt> GO TO **1 *4*553 

TOTFT=TOTFT*PLTFT * A*65* 

CALL 03PLT *4*655' 

**l CONTINUE *4*656 

TSMPaTSMP*DELP a 557 

*50 CONTINUE a 553 

C A 65 ? 

C INCREMENT TIME TO NEXT STOPPING POINT. ft 660 

6 IF TIME IS GREATER THAN TmX - STOP CALCULATIONS. ft 661 

C IF OIL FRONT HftS BROKEN THROUGH IN AT LEAST ONE 'STREAMLINE ft 662. 

C *NCOUNT c i .SET DELT * DELT*10. ft 663 

C A 66*. 

IF(NCOUNT.EQ.l) DELT=OElTT A 665' 

TSTOP«TSTOP*DELT A 666. 

IF (TST0P.GT.TMX) GO TO *60 A 667 

GO TO 230 A 668 

*60 CONTINUE A 66? 

CALL ENDPLT *4*570 

*70 CONTINUE A 671 

*80 FORMAT <//. 8 X»*NP * •» 13 ** NI * *.I3»* N 8 W * *«I3»* N 8 P e *.l 3 ft 57? 

1.* NST 3 * . 1 3 ) A 573 

*90 FORMAT (//.8X.*H « *.Fl0.5,* DElT 3 **Fl0.5»* TMX » **Flo.l** A 67* 

1RI 3 *.F10,5) A 575' 

500 FORMAT </7.8X,*DELP * *»Fl0»5t* *MX 3 *,Flo t l.* Y MX * *»Flo.l» A 576 

1 * QMNP * * . FI 0 , 5 ) A 677 

510 FORMAT <//.8X.*POR * *.F6.5»* SC* 3 *»F6.5'** SOI 3 *»F6.5»* So ft 573 

1R « *,F6,5.* SGR ■ *.F6.5) A 67? 

520 FORMAT <//*8X,*0ILl MOB 3 *,FiO,5,*' WAT MOB* *»F10.5»* GAS M03' 4 69o 

l« *.F) 0.5) A 631 
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530 FORMAT (6F10.0) A 092 

540 FORMAT (6110) A 693 

550 FORMAT (6F10.2) A 594 

560 FORMAT ( 2F1 0 • 0 > A 595 

570 FORMAT (*NO OIL BANK WILL FORM BECAUSE THE INITIAL! OIL SATURATION A 595 

1 IS NOT GREATER THAN THE RESIDUAL OIL SATURATION.*) A 697 

580 FORMAT (12F5.0> ' A 599 

590 FORMAT C 1 Ml ) ' A 69? 

595 FORMAT <8X. •INTERNAL WELL BALANCE DEFEATED*) A 599 

600 FORMAT ( 8X t "BOUNDING COMPLETE*) A 591 

610 FORMAT ( 8X * *RUN9ACK COMPLETE*) A 69>, 

620 FORMAT <//»SX»*WELL. NO *,I2,* AT **E10.3»* ♦ • * El 0 . 3 . *.*) A 593 

630 FORMAT (/.10*»*GAS RATE a *,Elo.3**. OIL RATE * *,FB,2**. WATER R A 594 

1 ATE = **F8.2**.*) A 595' 

640 FORMAT </*10Xt*CUM OIL * *»E10.3,*. CUM WAT » *»E10.3**. GOR a A 595 

1 E10.3»* . WOR * **E10.3) A 697 

650 FORMAT (/15X,«FIEL0i OIL RATE = *,£10.3,*. WATER RATE * **E1o.3) A 599 

651 FORMAT </15X,«FIEL0 CUM OIL a *,El0.3*** FIELO CUM WATER a *, A 59? 

1E10.3) A TOO 

660 FORMAT (/15X,«TIME a *»ElO*3.*DAYS. CONDUCTIVITY RATIO * **E10.3 A 70 1 
1 *.*.///) A 702 

END A 7 o 3 

8 1 

SUBROUTINE to BOUND IRREGULAR RESERVOIRS B 2 

REFl J.-K. LIN, PH.D. DISSERTATION* THE UNIVERSITY OF 83 

TEXAS AT AUSTIN (OECEMBER, )97 2). - 8 4 

8 5 

SUBROUTINE BOUND 8 f, 

COMMON/AAA/ XW (125) ,YW(125) *HI (125» *0(125) ♦ X3(125) » 8 7 

IY0(125)»NbP.NP,NI,N8W,H 8 9 

COMMON/BBb/ C< 55,55). 0(55) *OI<55) 8 9 

DIMENSION XI(55).YI(55)*XSB(100)»ySB(100)*XI3(100) , 8 10 

1Y!B(100>«AA(100*5S),88(100>*DS(100> 8 11 

N8P1=N8P*1 8 12 

PAI*3, 14159 8 13 

L«NP*NI 8 14 

00 10 J 5 l ♦ NSW 8 15 

XKJ)cXW(J.L) 8 16 

YI (J)*YW(J*LJ 8 17 

io continue a 13 

XBIN8PJ ) *XB ( 1 ) 8 1? 

YBIN8P1 ) sYB( 1 ) 8 20 

CC*0.5 8 21 

DO 20 1*1 , MSP 8 22 

XH«0.5*(X8<I*1)*XB(I>) 8 23 

VM*0.5*IVB(I»1)*Y8CI)) 8 24 

DSm=S0RT{(X8(I*l)-X8<I))**2*irBCl*l)-YB(I))**2) 8 25 

HDS*CC*USm 8 26 

TN»YBCI*1)-Y3(I) 8 27 

TD*X0(I*1)-X3(I) 8 29 

THETAaATAN2(TN,T0) 8 2? 

BEtA»THETA-0,5*PAI 8 3j 

XI8(I)aXM*H0S*C0S(BETA) 0 3l 

XS8II)aXM*H0S*C0S(BETA*P4I) 8 32 

YI8<I)*YK*H0S*SIN(BETA> 8 33 

YSB(I)sYM»H0S*SIN<6ETA«PAI> 8 34 

20 CONTINUE 8 35' 

NsNbv 8 36 

M«NBP 8 37 

00 30 1*1 *N8P 8 39 

DO 30 J*I»N8W 8 39 

FNUH»{XSB(I)-XIIJ))**2*(rs8(I)-YI(U))**2 8 *0 
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FDEN*<XIB( I) -XI ( J) ) »*2* < r I a ( I) -YI ( J) ) **2 B 41 

FACT=FNUM/FDEN B 42 

AA(1*J)»AL0G(FACT) B A3 

30 CONTINUE B 44- 

DO 50 J- 1 » N9P 8 45' 

$UM=0.0 B 45 

DO 4o K=1*L B 47 

FN»MXl3(J)-Xrf<K) )**2*(YIB(J)-YW(K) ) **2 B 43 

FDe(XSB(J)-XW(K) )*»2*(YSB<J)-t'W(K) )«»2 B 4? 

DSUMsQ (K ) /HI (K) ** ALOG (FN/FD ) B 53 

SUM=SUM*DSUM B 51 

*0 CONTINUE B 52 

BB ( J) cSUM 3 53 

50 CONTINUE B 54 

DO 7 0 I = 1»N B 55' 

DO 7 0 X=1*I B 55 

C<K»I)«0.0 B 57 

DO 60 J=1*M B 53 

C<K*I)=C(K*I)*AA(J,I)*AA(J»K) B 5? 

60 CONTINUE B 53 

C(I*K)=C(K*I) B 51 

70 CONTINUE B 52 

DO B0 I = 1*N B 53 

0 ( I ) = 0 « 0 B 54 

DO B0 J= 1 * M B 55' 

D(I)sD(I)*BB(J)*AA(J,I) B 55 

60 CONTINUE B 57 

CALL SOLVE (N) 8 53 

DO 90 J=1.NB# B 5? 

0(J*L)sOI«J)*H B 7 3 

90 CONTINUE B 71 

RETURN 8 72 

END B 73 

C 1 

SUBROUTINE solve USED By SUBROUTINE bound . SEE. BOUND C 2 

FOR REFERENCE. C 3 

C V 

SUBROUTINE SOLVE (N) C 5 

COMMON/BBB/ C<55»55) »0(55) iOI (55) C 5 

DIMENSION A ( 55 * 55 ) * ID ( 55 ) C 7 

NN=N*1 C 9 

MsN C 9 

DO 10 I = 1*M C 13 

A ( I *NN) =0(1) C 11 

ID ( I ) s I C 12 

DO 10 J«1»N C 13 

A(I*J)=C(I»J) Cl 4 

10 CONTINUE C 15 

K»1 C 15 

20 CONTINUE C 17 

NRsK C 13 

NC«N C 1? 

BoABS< A (K»K) ) C 23 

DO 40 I=K»M C 21 

DO 40 J=K*N C ?? 

IF (ABS(A ( I, J) ) -8) 40*40.30 C 23 

30 CONTINUE C 24 

NR«I C 25' 

NC«J C 25 

B»ABS( A(I.J) ) C 27 

40 CONTINUE C 23 

IF (NR-K) 70 * 70*50 C 2? 
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50 


continue 


c 


3D 




DO 60 J a K* NN 


w 


31 




CaA (NR. J) 


c 


3? 




A ( NR, J) =A ( K» J) 


c 


33 




A (K, J) sC 


c 


3<> 


60 


continue 


c 


35 


70 


continue 


c 


35 




IF (Nc-K) loo, 100.80 


c 


37 


60 


continue 


c 


39 




00 90 I=1*M 


c 


3? 




C=A ( I ,NC) 


c- 


4D 




A ( I » NC ) 3 A ( I » K ) 


c 


41 




A (I ,K) sC 


c 


42 


90 


CONTINUE 


z 


43 




Is ID ( NC ) 


c 


4 ^ 




ID(NC) slD(K) 


c 


4 5 




ID ( K ) = I 


c 


4 S 


100 


CONTINUE 


r* 


47 


110 


continue 


C 


49 




IF ( A ( K « K ) ) 120*180*120 




45 


120 


CONTINUE 


c 


5 D 




KK S K* 1 


r 

w 


51 




00 lAo JsNK.NN 


C 


5? 




A (K* J) sA (K» J) /A(K*<) 


C 


53 




DO lAO I=1*M 


c 


5 '4 




IF (K-I) 130* 1A0.130 




55 


130 


CONTINUE 




55 




A(I, J) sA(I, J)-A(I,K)*A(K,J) 


r* 


57 


140 


continue 


* 


59 




fc=KK 




59 




IF (K-N) 20*110*150 




5D 


150 


continje 


z 


51 




DO 17o I=1*N 


r 


62 




DO 170 Us 1 , N 


c 


S3 




IF (IO(J)-I) 170.160*170 


c 


6 '♦ 


160 


CONTINJE 


r 


6 5 




01 (l)sA(J.NN) 


c 


65 


170 


continue 


c 


57 




GO TO 190 


c 


63 


160 


continue 


r* 

s* 


69 




PRINT 200 


P* 

w 


7D 


190 


continue 


c 


71 




return 


c 


72 


200 


format (19H NO UNIUUE SOLUTION ) 


c 


73 




ENO 


c 


7* 



*P* l 

THIS SUBROUTINE plots the locations OF THE OIL and *P* 2 

water fronts as they dEvElope. all data is taken prom *p° 3 

THE MAIN PROGRAM. *P* A 

*P* 5' 

subroutine obplt *p* s 

COMMON/AAA/ XW (125) ,YW<125) .Hi < 125) «Q<125) , X8< 125) . *p* 7 

1YB1125) *NBP»NP,NI» NBW.H *P* 9 

COMMON/CCC/ X(80.20) .Y<80»2;)) .Xo<80.20) ,YO<90»20) * *P* 9 

1WQ(B0,20) *SCL*XMX»YMX»NST»OmNP,RGO<80*20> « *P* ID 

2SG0R.N9lw.TMEI *P* 11 

CALL SYMBOL < 0 . 5 , 1 0 . 5. . 1 A, 7hT IM£ = .27o,0.7) *P tt 12 

CALL NUMBER (999,0* 999.0 » »lA.TM£ 1, 270.0.0) *p* 13 

CALL SYMBOL (999,0.999.0*. U.5H OA YS * 270 . 0 . 5 ) *P* 1A 

CALL PLT (1.125*XMX*SCL*l.l25*YMX*SCL*-3) *P* 15 

CALL PIT (XBlNBPJwSCLtYbtNdPlwSCL.B) *P* 15 

00 10 1 3 1 * NBP *P* 17 
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CALL PLT (X8(I)<*SCL.YB(I)*SCL*2> «P* 13 

10 CONTINUE *P* 13 

00 30 I a 1 * NP I W 20 

IF (0(1) .GE. 0.0) GO TO 20 *P* 21 

CALL SYM80L'(XW(I)*SCLiYK(I)*SCLt.07tl.o*-l) «P* 22 

GO TO 30 *P* 23 

20 CONTINUE «P* 24 

CALL SYM80L (XW < I ) *SCL • YW < I > *SCL * . 07 * 7 * 0 • -1 ) *P* 2 5 

30 CONTINUE *P« 23 

00 50 J = 1 * NP I W *P* 27 

IF (Q(U) .oE.0.0) GO TO 5o *P* 23 

NSL a -NST*Q (0) /QMNP*0.5 *P* 23 

00 50 <=1*NSl a P* 30 

IF <A8S(X0(J»K) ) .GT.XMX) 30 TO 40 *P* 31 

IF (A8S ( YO ( J*K) ) .GT.YMX) GO TO 40 *P* 32 

IF (RGO(JiK) .LE.SGOR) GO TO 40 *P* 33 

CALL SYMBOL (XO( Ji<)*SCLiYO(J,K)*sCLt.07»31»0»-l) *P* 34 

40 CONTINUE *P* 35 

IF (A8S(X(J.K) ) .GT.XMX) GO TO 50 *p* 35 

IF (A0S(Y (J.K) ) .GT.YMX) GO TO 50 *P* 37 

IF (WQ( J,K) ,GT. 0.0) GO To 50 *P* 33 

CALL SYMBOL <X(J.K)*SCL»Y(J.<)«SCL».07.40t0t-l) °P* 33 

50 CONTINUE <*P* 40 

CALL PLT (0.0*0.0*999) *P* 41 

RETURN " «P a 4? 

ENO *P« 43 
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